The properties of amorphous GaN 
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Abstract: In this paper, we present three amorphous GaN models obtained from the first princi- 
ples simulation. We find that a chemically ordered continuous random network is the ideal structure 
for a-GaN. If we exclude the tail states, we predict a 3.0eV optical gap for 64-atom model and 2.3eV 
for 250-atom models. We observe a highly localized valence tail and a remarkably delocalized expo- 
nential conduction tail which we associate with different hybridization in the two tails. Based upon 
these results, we speculate on potential differences in n and p type doping. The structural origin 
of tail and defect states is discussed. The vibrational density of states and dielectric function are 
computed, and are consistent with experiment. 



With a profound impact on lighting technology and 
other applications, crystalline GaN bas been the sub- 
ject of vast inqunyi - — . However, lattice mismatch with 
substrates makes it difficult to grow. Recently, amor- 
phous GaN has become attractive due to its potential to 
solve the lattice match problem and its natural isotropy. 
A number of experiments have investigated the struc- 
tural and optoelectronic properties of a-GaN 4 9 . In 1997, 
we proposed that a-GaN might find use as an electronic 
materia l 10 ! 11 . 

The structure of a-GaN is controversial. Some ex- 
periments observe a large concentration of homopolar 
bondsi^, contradicting other studies 8 . We did not ob- 
serve N-N or Ga-Ga bonds in our computer models, but 
other calculations suggest a more disordered network^. 
Doping of a-GaN is an important topic, but a full under- 
standing of the intrinsic electronic features of undoped 
a-GaN, like the origin of defect states and tail states, is 
a necessary precursor. Naturally, the detailed properties 
depend upon the mode of growth of the material; our 
work is most relevant to least defective "ideal" a-GaN. 

It has been more than a decade since the first a-GaN 
model was proposed, a time during which many experi- 
ments have been carried out, and simulation tools have 
experienced major developments. Therefore, additional 
calculations have been undertaken to generate atomistic 
models of a-GaN and to further explore its interesting 
properties. 

In this paper, we propose atomistic a-GaN models 
formed via ab initio molecular dynamic simulation with 
a plane wave basis. The network topology is analyzed 
through radial and angular distribution functions, and 
structural statistics. We find that Ga and N atoms 
strongly prefer to be four-fold, and homopolar bonds are 
rare in the network. We also predict electronic proper- 
ties and connect the electronic structure to the topology 
of the network. We show that the conduction edge has 
Urbach (exponential) form and is extraordinarily delo- 
calized, and the valence edge is very sharp with highly 
localized states. Doping is briefly discussed. Vibrational 
properties and dielectric functions are predicted. 

All calculations in this work are performed with the 
Vienna Ab-initio Simulation Package (VASP)^ based on 
density functional theory (DFT) within the local density 



approximation (LDA) and Vanderbilt's ultrasoft pseu- 
dopotentials. Three models are generated via computer 
alchemy: a 64-atom model (obtained from a 64-atom InN 
model of Ref. 15), 250-atom model I (obtained from 250- 
atom a-GaAs model of Ref. 1^) and 250-atom model 
II (obtained from 250-atom a-InN model of Ref. [TBI ). 
The 64-atom model is annealed and equilibrated at 500K 
and 250-atom model I is annealed and equilibrated at 
300K. After zero-pressure relaxation, the density of the 
64-atom model is 5.8g/cm 3 , 95% of the crystalline GaN 
density; both 250-atom models have mass density around 
5.6g/cra 3 , 92% of the crystalline GaN density. The new 
64-atom model shows a improved cohesive energy (lower 
by 0.16eV/atom) compared with the previous model 10 . 
Where chemical order is concerned, there are no homopo- 
lar bonds in the 64-atom model and 250-atom model I, 
and only one N-N bond in the 250-atom model Has one 
might expect for a partly ionic system. For coordination, 
we note that most atoms tend to be four-fold, suggest- 
ing that a-GaN retains the zinc-blende/wurtzite char- 
acter of crystalline GaN, in significant contrast to our 
early study 10 . We list the structural properties in Table 
H Since the 64-atom a-GaN was made from melt-and- 
quench method, while the both 250-atom models are gen- 
erated only by a energetic relaxation, the little increase 
of over- and under-coordinated atoms in both 250-atom 
models could be the artifact of relaxation technique used. 
This result suggests that an ideal a-GaN is a chemically 
ordered continuous random network. 

We plot the partial pair-correlation functions of the 
three models in Fig. [TJ For the Ga-N partial, there exists 
a sharp first peak around 1.94A for all three models. The 
pre-peak/shoulders in Ga-Ga partial (marked by black 
arrow) indicate that there exists two local environment 
for Ga atoms and we will show that those two Ga sites are 
related to edge-sharing N tetrahedral structures. Due to 
the N-N bond, there is a small peak around 1.53 A in N- 
N partial of 250-atom model II (marked by black arrow). 
Overall, the pair correlation functions of three models 
exhibit similar features and they are close to the results 
of Ref. HI 

Next, we analyze the angle distributions for Ga-N-Ga 
bonds and N-Ga-N bonds. The 250-atom models yield 
major peak positions close to Ot = 109.47° for both Ga- 
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TABLE I. The statistical distribution of the main structural components of three models. 
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FIG. 1. (Color online) Partial pair correlation functions of 
the three models (see text). In the Ga-Ga partial, the pre- 
peak/shoulders are marked by a black arrow. For the N-N 
partial, the black arrow marks the peak due to the only N-N 
bond in 250-atom model II. 



N-Ga and N-Ga-N angle. For the 64-atom model, the 
major peak positions are slightly off #t, the mean value 
of of N-Ga-N and Ga-N-Ga angle, being 109.15° and 
108.65° respectively. Thus, we conclude that a-GaN re- 
tains strong vestiges of its crystalline short-range order 
and tends to form a tetrahedral structure. We observe 
a pre-peak around 80° for the 64-atom model (shoulders 
for the 250-atom model) in the Ga-N-Ga angle distribu- 
tion which implies that there are two distinct sites for Ga 
atoms. After a detailed investigation, we find that the 
small angle is due to edge-sharing units with distorted 
angles (appearing as four-member-rings with Ga-N-Ga 
angle between 75° and 95°). We will show that this kind 
of distortion is responsible for some of the electronic tail 
states. 

The vibrational properties of a-GaN are characterized 
through the vibrational density of states (VDOS). Start- 
ing with a thoroughly relaxed cell, the force constant and 
dynamical matrix is obtained from perturbing each atom 
in turn by 0.015A, and computing forces on all atoms in 
the model for each perturbed conformation. The VDOS 
of the 64-atom model and 250-atom model II are reported 
in Fig. [21 Both models show similar features. For com- 
parison, we also plot the VODS of crystalline GaN from 
Ref. LL7| as insert. Our results show that the amorphous 



FIG. 2. (Color online) Vibrational density of states of 64- 
atom model and 250-atom model II. The eigenvalues were 
Gaussian broadened with a width of lmeV. The VODS of 
crystal wurtzite GaN is plotted as an inset from Ref. 



VDOS retains some features of crystal VDOS such as the 
two peaks in the first band. However, we did not observe 
two distinguished peaks in the optical band 18 , and the 
gap between the acoustic band and the optical band fills 
in substantially. The results are quite consistent with a 
recent Raman study 8 . 

We describe the electronic structure by analyzing the 
electronic density of states (EDOS), inverse participation 
ratio (IPR) of the individual states, and dielectric func- 
tions. Figj3] shows the EDOS of the three models. Over- 
all, all three EDOS have similar character with slight 
differences in detail. The conduction band tail is dra- 
matically broader than the valence band tail. If we de- 
fine the band gap as the difference between highest ex- 
tended energy level of states in the valence band and 
lowest extended energy level of states in the conduction 
band (excluding the mid-gap and tail states), we esti- 
mate that the band gap to be about 3.0eV for the 64- 
atom model and 2.3eV for the 250-atom models. We 
indicate this in Fig|4] (the region between black arrows 
is taken as optical band gap). The band gaps obtained 
from our models are smaller than the experimental value 
3.1eV— . This is primarily due to the LDA which is well- 
known to underestimate the band gap 19 . By fitting the 
conduction band tail to an exponential, we report the Ur- 
bach energy, E u « 420raeV for the 64-atom model and 
E u w 490meV for the 250-atom model I, comparable to 
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FIG. 3. (Color online) Electronic density of states of 64-atom 
model(32 K-points are used), 250-atom model 1(8 K-points 
are used) and 250-atom model 11(8 K-points are used). The 
Fermi level is at eV. 



1 1 
1,1 


a 


i ■ i ■ i 1 i ■ 

b 

1 1 64-atom 

i 

, U , ■ 1 1 1 ■ 1 


-1 

' iiJlillll 


c 


d 1 2 3 

250-atom model 1 

r . . ) 


-1 g 

; i 

.Hilhll, 


■ ■■■■■ 
f 1 2 

250-atom model 11 

.1 ■■ ■ 



1 

Energy(eV) 



FIG. 4. (Color online) Inverse participation ratio analysis for 
three models. Large IPR implies strong localization. For each 
model, the optical band gap are estimated by excluding the 
mid-gap and tail states (region between two black arrows in 
each plot). Note the highly delocalized conduction tail. The 
Fermi level is at OeV. 

the reported value "several hundred meV" in Ref. [20|. In 
addition, for 250-atom model II, there are defect states in 
the deep band region between -18eV and -12eV far below 
the Fermi level. These defects are due primarily to the 
N-N bond. Other than the defect states in deep band 
below -8eV, we did not observe any electronic signature 
of N-N bond around the optical gap. 

To characterize the localization of the tail states 
around the gap, we performed an IPR analysis for all 
three models. The IPR measures the degree of local- 



ization given an electronic state 21 . For highly localized 
states, IPR=1; for extended states, IPR=1/7V, where N 
is the number of atoms. The results are plotted in FigHJ 
By projecting the EDOS onto different atomic orbitals, 
we find that the valence tail is built from N-p, Ga-p and 
Ga-d orbitals. This implies a high sensitivity to bond 
angle disorder, which is presumably the reason for high 
localization. The conduction tail is localized on Ga-s 
and N-s orbitals. Since the s-s interaction is only affected 
by bondlength, the conduction tail states exhibit remark- 
ably weak localization and the conduction tail is almost 
immune to angle disorder. This situation is somewhat 
similar to a-Si02, where there is also a large asymmetry 
on IPR between valence band and conduction band tails 
as discussed by Robertso n 22 ! 23 . To our knowledge, this 
effect has not been reported in nitrides. The asymmetry 
in width and localization of tail states suggests that n 
and p doping for a-GaN will be quite different 24 . Due 
to the highly localized valence-band tail states, it will be 
more difficult to move the Fermi level toward the valence 
mobility edge, complex compensations may happen and 
mobility is likely to be poor. Thus, the p-type doping 
is expected to be relatively more difficult than n-type 
doping to obtain the same carrier concentration 2 ^. 

To correlate electronic structure with topological units, 
we picked seven electronic states (a-g in FigHJ with rela- 
tively high IPR and projected them onto individual atom 
sites. In FigEJ we present the characteristic atomic sites 
associated with those tail states a-g. In the 64-atom 
model, the tail states a and b are highly localized on 
atom- 13 (FigJSfa)) whose four neighbors are almost in 
the same plane, with distorted Ga-N-Ga angle 75°; the 
state b is also localized on atom-28(Fig|5jb)), the only 
three-fold N atom in the network which formed a small 
Ga-N-Ga angle near 89°. In 250-atom model I, three- fold 
N atom-43 (FigGJc)), which formed a 87° Ga-N-Ga an- 
gle, is strongly associated with eigenstate c and d; the 
tail state d are also localized on atom- 70 (Figj5jd)), a 
three-fold N atom with all its neighbors almost in the 
same plane; moreover, the conduction-band tail state e is 
mainly localized on atom- 162 (Figj5fe)), a three- fold Ga 
atom with disordered N-Ga-N angle. In 250-atom model 
II, three- fold N atom-Ill whose three neighbors are al- 
most in the same plane (FigJSJf)), contribute more to the 
valance-band tail state /; two four- fold N atoms, atom-4 
and atom- 50, formed edge-shared tetrahedron with disor- 
dered Ga-N-Ga angle, are strongly associated with elec- 
tronic state g. Overall, atoms with distorted angle are 
associated with valence tail states. Finally, we briefly re- 
mark that, for the 64-atom model, the imaginary part of 
dielectric function e(w) has a major peak position around 
6.8eV for all three directions. This result is comparable 
to the experimental work reported in Ref. 

In conclusion, we created a-GaN atomistic models us- 
ing state of the art methods. Most atoms in the network 
tend to be four- fold and form tetrahedral structures. We 
predict a 3. OeV optical band gaps for 64-atom model and 
2.3eV for 250-atom model. We find an interesting and 
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FIG. 5. (Color online) Atomistic origin of electronic tail and gap states correlated with state a-g indicated in FigUJ Dark(small) 
atom is N; light (large) atom is Ga. (a) atom- 13 is associated with tail states a and b. (b) atom-28 is associated with tail state 
b. (c) atom-43 is associated with tail states c and d. (d) atom- 70 is associated with tail state d. (e) atom- 162 is associated 
with tail state e. (f) atom-Ill is associated with tail state /. (g) atom-4 and atom-50 are associated with tail state g. 



large asymmetry in localization between valence and con- 
duction tail due to the different orbital interaction, which 
should yield quite distinct properties in n and p type 
doping. The atomistic origin of tail and defect states is 
discussed, and the disorder in bond angle is likely to in- 
troduce valence tail states, whereas the conduction tail is 
due primarily to bond length disorder. The vibrational 
density of states retains some qualitative features from 
the crystal, and the dielectric functions shows a peak 



around 6.8eV, both of which are in agreement with ex- 
periment. Our work focuses primarily on "ideal" GaN to 
establish a reference model. Ion bombarded samples are 
indeed likely to exhibit far more disorde r 12 i 13 . 
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